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ABSTRACT 

We describe a new program called cryptic splice 
finder (CSF) that can reliably identify cryptic splice 
sites (ess), so providing a useful tool to help inves- 
tigate splicing mutations in genetic disease. We 
report that many ess are not entirely dormant and 
are often already active at low levels in normal 
genes prior to their enhancement in genetic disease. 
We also report a fascinating correlation between the 
positions of ess and introns, whereby ess within the 
exons of one species frequently match the exact 
position of introns in equivalent genes from another 
species. These results strongly indicate that many 
introns were inserted into ess during evolution and 
they also imply that the splicing information that lies 
outside some introns can be independently recogn- 
ized by the splicing machinery and was in place prior 
to intron insertion. This indicates that non-intronic 
splicing information had a key role in shaping the 
split structure of eukaryote genes. 

INTRODUCTION 

Eukaryotic genomes contain large numbers of splice sites, 
known as cryptic splice sites (ess), which are generally held 
to be disadvantageous sites that are dormant or used only 
at low levels unless activated by mutation of nearby au- 
thentic or advantageous splice sites (1,2). Once activated, 
ess may be used extremely efficiently, resulting in a wide 
range of genetic disease (3-5). It is generally accepted that 
ess are suppressed by nearby stronger splice sites and that 
splice site selection can be viewed as a competition 
between the various potential splice sites in a pre-mRNA 
for the splicing machinery (1,2). 

For genes with many introns it is suspected that up to 
50% of mutations that cause disease do so by affecting 
splicing, either through the activation of ess, exon 



skipping or disruption of alternative splicing (4-7). Css 
are found in exons as well as introns and their recognition 
by the splicing machinery is similar to splice site recogni- 
tion in general and is dependent upon information both at 
the splice site and outside this region at enhancer and 
silencer sequences (8-10). 

It is important to be able to predict the positions of css 
that might be activated in genetic disease and a number of 
DNA-sequence scanning programs have been developed 
for this purpose. Such programs are often highly inform- 
ative but are handicapped by the complex nature of the 
nucleotide information that is required to define a splice 
site (4,8,11,12). 

Our previous work indicates a connection between css 
and introns. We identified a small number of css in the 
exon regions of actin genes by experiment and discovered 
that eight out of nine of these exonic css sites exactly 
match the positions of introns in actin genes from other 
species, which led us to conclude that these particular 
actin introns were inserted into css during evolution 
(13,14). This finding may help to explain why and how 
eukaryotes acquired introns; however, it is important to 
establish if our results for the actin gene family are gener- 
ally applicable. 

We have been unable to identify a DNA-scanning 
program that reveals the same strong correlation between 
predicted actin css and intron positions that we observed 
through experiments (13). However, this is probably 
because DNA-scanning programs were not designed spe- 
cifically for this purpose and because of the difficulties 
such programs face in distinguishing between css and 
false-positive non-functional splice sites (12). Here, we de- 
scribe a program called cryptic splice finder (CSF) that can 
reliably identify css by EST-to-genomic alignment. It does 
this by identifying transcripts that have been generated 
through the low level use of css by normal genes. Unlike 
the scanning programs, CSF cannot predict the position 
of splice sites that are created de novo by gene mutation. 
However, this program provides a useful complementary 
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resource for studies of genetic disease and it also enabled 
us to establish that there is a strong and general correl- 
ation between the positioning of ess and of introns. The 
evolutionary implications of this finding are discussed. 



MATERIALS AND METHODS 

CSF 

CSF (http: //www. ncbi.nlm.nih.gov/IEB/Research/csf) 
predicts ess based on spliced alignment of ESTs. Each 
EST is aligned against the genome independently from 
other ESTs. In order to be considered by CSF, a gap in 
the alignment must be flanked by a minimum number (25 
by default) of matching residues. CSF searches for EST 
alignments that form the patterns illustrated in Figure 1 A. 
It can be seen that the majority of ESTs must share a 
common gap or deletion and in addition must include 
minor transcripts that share only one of the common 
deletion endpoints. CSF defines the common deletion 
endpoints as authentic splice sites and the deletion end 
point of the minor transcript(s) as cryptic or alternative 
splice sites (arrowed). For constitutively spliced genes, 
these minor deletion endpoints occur very infrequently 
and are therefore candidate ess. For alternatively spliced 
genes, CSF identifies both ess and more frequently 
occurring alternative splice sites (see 'Results' section). 
As illustrated in Figure 1A, ess can be 5' or 3' and 
upstream or downstream of the alternative authentic 
splice site. In the majority of cases the authentic splice 
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Figure 1. (A) CSF searches for transcript alignments that form one of 
four patterns (i-iv). All of these patterns contain a group of major 
transcripts that share a common deletion and a minor transcript that 
shares only one of the deletion endpoints. CSF defines the common 
deletion endpoints as authentic splice sites and the less common 
deletion endpoint of the minor transcript(s) as cryptic or alternative 
splice sites (arrowed). (B) Schematic of the HBB gene for human 
P-globin, which contains two introns that are constitutively spliced 
from pre-mRNA. As illustrated the vast majority of ESTs align as 
shown and define the three exons of this gene. The circle shows a 
pattern of ESTs that CSF is designed to recognize and that is 
reported in Figure 2A. The numbers in brackets show the genome co- 
ordinates of the three splice sites that are identified and listed by CSF 
(Figure 2A). 



sites defined by CSF are the same as the splice sites 
defined by reference sequences (i.e. NM_000518.4) and 
so represent commonly used splice sites (see below for 
exceptions). 

In more detail, coordinates of the splice sites from 
adjacent exon-intron-exon sequences are pooled into 
four-tuples which are then loaded into a relational 
database alongside the data linking them to their align- 
ments. The database runs a query to detect tuple pairs 
satisfying the ess condition: the overlapping introns 
match at one end and mismatch at the other, with a mis- 
matching intron end from one tuple residing within the 
exon from the other tuple. For each tuple pair returned 
by the query, the splice sites of the intron that has more 
supporting ESTs are declared authentic and the remaining 
site is declared cryptic. Splice site coordinates are mapped 
against the NCBI36/hgl8 human genome assembly. 

As illustrated (Figure 1A), a ess detected by CSF is fur- 
ther classified as to whether it is a 5' or 3' ess and whether 
it is located in an exon or intron. The number of tran- 
scripts supporting a splice (either authentic or ess), is 
printed and the complete list of such ESTs is linked 
under the count. Where applicable, the count is followed 
by a number in parentheses. For exonic ess, parentheses 
always appear on the left-hand (authentic) side of the CSF 
report and on the right-hand side (cryptic) for intronic ess. 
The numbers in parentheses show the transcripts that 
formally satisfy the ess conditions (see Notes about CSF 
section). 

Notes about CSF. CSF only lists authentic 5'- and 
3'-splice sites when one of the two authentic splice sites 
is a more common alternative for another nearby splice 
site that is listed under the cryptic (alternative) column. 
This means that CSF only lists a subset of the total 
number of authentic splice sites for any one gene. 

CSF provides a very good way of screening large 
numbers of transcripts for the minority that are likely to 
have been generated by use of a ess or nearby alternative 
splice site. However, candidate transcripts do need to be 
checked, by for example, using the Splign alignment option 
(15). CSF classifies a minority of less well-supported dele- 
tions as authentic splices because these happen to meet 
the CSF conditions. These sites can be easily recognized 
from the CSF output because they have lower levels of 
support than other authentic splice sites and do not 
match reference sequence splice sites. Caution needs to 
be taken with the interpretation of these particular 
'authentic' splice sites and also with the interpretation 
of ess that are paired with authentic splice sites that 
have a low-parenthesis score. See Supplementary Data 
for illustrative examples. 

Searching CSF. Two methods of searching for ess are 
provided (Supplementary Figure SI). In the first method, 
a landmark EST accession is submitted and CSF returns a 
list of ess within the genomic range of the entered EST. 
In the second method, an arbitrary genomic interval can 
be specified and CSF will return a hierarchic list of ess for 
that interval. Entire chromosomes can also be entered 
such as NC_000001 (human chromosome 1) for which 
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CSF currently lists 3232 ess. Links are provided in order 
to view the sequence alignments of individual ess. CSF can 
predict ess for Homo sapiens, Bos taurus, Mus mus cuius, 
Danio rerio and Arabidopsis thaliana and we intend to 
expand this range as further transcript data becomes 
available 

Statistical analysis. There are 91 known different intron 
positions within the coding region of the actin gene family, 
these have been identified by sequencing actin genes 
from over 160 different species (16-18). The positions of 
all of these introns together with the 14 ess that we have 
identified are plotted in Supplementary Figure S2A. Actin 
genes usually have 375 codons and therefore three times 
this number of possible positions for introns or ess and so 
the probability of a single ess exactly matching an intron 
position by chance is 91/(3 x 375). The exact Fisher's test 
gives the P- value of 1.6x 10 -10 for 11 or more matches 
out of 14 occurring by chance. Similarly, we identified 135 
ess in the coding region of 51 different genes of the ribo- 
somal protein gene database. The 51 genes are 190 codons 
in size on average and have a total of 957 introns (from 
22 species). The probability of a ess exactly matching an 
intron by chance is therefore 957/(190 x 3 x 51) = 0.033. 
The binomial probability of 33 or more matches out of 
135 occurring by chance is P = 2.2 x 10 -15 . 

RT-PCR. Total RNA was extracted using Trizol and the 
cDNA first strand was synthesized from 0.1 jig of RNA 
using superscript III (Invitrogen) and random hexamers. 
PCR products were generated with Taq polymerase (NBI) 
for 25 or 35 cycles and separated on 5% native polyacryl- 
amide gels. PCR bands were excised and cloned into 
pGEM-T Easy vectors (Promega) for sequencing by 
colony PCR followed by ABI Prism Big Dye Terminator 
cycle sequencing (Applied Biosystems). 

RESULTS 

Principle of ess detection by CSF 

Css are often used highly efficiently in genetic disease fol- 
lowing the mutation of nearby more competitive splice 
sites. We, therefore, reasoned that css might be used at a 
low but detectable frequency by normal genes. Figure 1A 
shows the patterns of EST alignments that CSF is designed 
to identify. It can be seen that CSF identifies groups of 
ESTs or transcripts that share a common deletion, when 
aligned to the genome, together with a minor transcript(s) 
that shares just one of the common deletion endpoints. 
CSF defines the common deletion endpoints as authentic 
splice sites and the unusual deletion endpoint of the minor 
transcript as a cryptic or alternative splice site (see 
'Materials and Methods' section for further details). 
Figure IB illustrates how ESTs align to the human HBB 
gene for (3-globin, which has two introns. The circle 
identifies a pattern of alignments that is recognized by 
CSF because it includes a minor transcript that has an 
unusual deletion endpoint (position 5 204 752) that is a 
predicted css (Figures IB and 2A). The reason why such 
css predictions turn out to be accurate (see below) is 



because of the restriction that the predicted css is paired 
with a commonly used splice site. This distinguishes 
mRNA deletions that are generated by low-level aberrant 
splicing from deletions that are generated by non-splicing 
mechanisms such as errors during transcription or during 
the generation of the EST. 

Figure 2A shows the CSF output for HBB, which is one 
of the first genes in which css were identified (19). CSF 
identifies a single EST called BUI 98 526 as having been 
generated by aberrant splicing. BU198526 aligns to HBB 
as illustrated in Figure lB(ii), this alignment can be viewed 
in detail by a link to Splign (Figure 2). As explained 
above, the reason why BUI 98 526 was identified by CSF 
is because it forms a pattern alignment with other ESTs as 
illustrated in Figure lA(ii). The CSF output (Figure 2A) 
shows that BUI 98 526 shares a deletion endpoint at 
position 52044 605 that is in common with 703 other 
ESTs (link provided) but differs in having an unusual 5' 
deletion endpoint at position 5 204 752, which is the pre- 
dicted css. A comparison with the known css of HBB 
confirms that 5 204 752 is indeed a css (Supplementary 
Table SI). We show below that css predictions by CSF 
are very reliable even if supported by just a single EST, 
as in this case. However, HBB has 11 known css 
(Supplementary Table SI), which illustrates that CSF is 
limited in its predictions, most probably by the amount of 
available transcript data. 

For alternatively spliced genes, CSF identifies both css 
and also nearby alternative splice sites that act to shorten 
or lengthen exons. This is because the two types of splice 
sites are rather similar. For example, the css identified in 
HBB (Figure 2A) might be considered an alternative splice 
site if used at a greater frequency than 1 in 703. 

WT1 is an example of gene that encodes at least three 
alternatively spliced mRNAs and for which there are only 
80 ESTs. CSF identifies two css (Figure 2B), however, the 
CSF output also shows that the predicted css at 
32 370 103, for example, is supported by five ESTs 
including the reference sequence NM_000378.3 and that 
the 'authentic' splice site nine bases away at 32 370 094 is 
supported by only eight transcripts. The similar usage of 
the authentic and css shows that these are really alterna- 
tive 5' ss. However, it is useful to have this type of infor- 
mation because alternative splice sites are also implicated 
in genetic disease and in this particular case disruption of 
the splice site at position 32 370 094 gives rise to Frasier 
syndrome, possibly due to the increased use of splice site 
32 370 103 (20,21). 

To test CSF we analyzed a database called DBASS 
(database of aberrant splice sites), which lists 340 human 
genes that have one or more css that are activated in 
genetic disease (3). There are 814 different css listed in 
DBASS and CSF predicts 609 css from the same set of 
340 genes. Fifty-eight percent of these predictions are sup- 
ported by only single ESTs (Supplementary Table SI). 
Before comparing the css identified by CSF with those 
of DBASS, we first asked whether the css predictions 
that are supported by just single ESTs were likely to 
have been generated by aberrant splicing. We reasoned 
that because CSF identifies deletion endpoints irrespective 
of their sequence, then if these rare deletion endpoints 
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Figure 2. (A) CSF output for the human gene HBB for |3-globin. (B) CSF output for WT1 (see text). The coordinates that are used refer to the 
NCBI36/hgl8 human genome assembly. It should be noted that HBB and WT1 genes align in a 3'— 5' direction with respect to their genome 
coordinates. 



were generated by splicing they should look like splice 
sites. The predictions were divided into 5' and 3' ess and 
Table 1 shows that the predicted ess have a very good 
match to the expected 5'- or 3'-splice site consensus se- 
quence, which is a strong proof that the vast majority of 
even the rarest deletions identified by CSF were generated 
by aberrant splicing. 

Of the two sets of DBASS and CSF ess, only 46 are in 
common (Supplementary Tables SI and S2). However, 
there are only 61 cases where the DBASS and CSF ess 
are located within the same exon or intron 
(Supplementary Table SI, column 5), giving a match 
rate of 46/61 or 75%. Thirteen of the 15 DBASS ess 
that did not match were either created de novo by 
mutation and would, therefore, not be expected to be 
identified by CSF or were of opposite type, for example, 
a 5' ess in DBASS and a 3' ess in CSF and would, there- 
fore, be unlikely to match (Supplementary Table SI). Only 
two out of 48 ess predictions that could be meaningfully 
compared with DBASS did not exactly match 
(Supplementary Table SI). This result together with the 
clear consensus sequence results shown in Table 1 indi- 
cates that CSF predictions are highly reliable and can 
only generate a very low level of false positives even 
when supported by single ESTs. The small number of 46 
ess in common between the CSF database and DBASS 
(Supplementary Table SI) presumably indicates that 
both DBASS and CSF have identified relatively few of 
all possible ess within this set of 340 genes. 

We confirmed that ess predicted by CSF could also be 
detected by experiment. Five ess predictions by CSF that 
DBASS shows to be activated in patients were detected in 
human cell lines that do not have the causative genetic 
mutations (Figure 3, lanes 1-5). The last lane shows an 



example of a ess prediction by CSF that has not been 
reported in patients. The PCR products marked with as- 
terisks were confirmed by sequencing (Supplementary 
Figure S3). These results confirmed the CSF predictions 
and also show that the ess we analyzed were already active 
in a normal genetic background. 

Css were previously reported to be active at very low 
levels in normal globin genes (22) and our results system- 
atically extend this finding. It is clear from Figure 3 that 
some splice sites that have been classified as css are used 
by normal genes at a relatively high frequency (lanes 2-5), 
which is in accordance with the CSF analysis 
(Supplementary Table S2). These css can therefore also 
be regarded as alternative splice sites that are further 
activated in disease. 

Css and introns 

We previously identified nine css within the coding region 
of the actin gene family by experiment and reported an 
eight out of nine exact match to the position of introns in 
actin genes from other species (13,14). CSF identified eight 
css within the coding region of the actin gene family of 
which three were identical in position to the previously 
identified actin css. Of the five new css identified by 
CSF, three exactly match the position of introns in other 
species, therefore extending our previous study and adding 
further support to the validity of CSF predictions 
(Supplementary Figure S2). 

Because CSF predictions are sufficiently accurate, we 
could use this program to establish whether there is a 
general relationship between the position of css and 
introns. We therefore analyzed an extensive intron 
database compiled of some 80 genes from a wide range 
of species that encode ribosomal proteins (23). Figure 4A 
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shows a small part of this database that records the pos- 
itions of introns against a 14 amino acid stretch of the 
gene RPS5 from 23 species (file kindly provided by N. 
Kenmochi). Introns that fall between codons are indicated 
by red shading and introns that fall within codons are 
indicated in blue (Phase 1) and green (Phase 2). For 
example, there is an intron located within a codon for 
valine (V) in the fungus Cryptococcus neoformans but in 
no other listed species. This observation is typical in that 
although introns are wide-spread among eukaryotes, 
the vast majority of individual intron positions are 
found in only a minority of species (Supplementary 
Table S3) (18). 

We screened for ess in ribosomal protein genes from 
human, mouse, Danio rerio and Arabidopsis thaliana and 
identified 74, 41, 16 and 4 ess, respectively, within coding 
regions, these diminishing numbers simply reflect the 

Table 1. The alignment of ess predictions by CSF reveals consensus 
sequences typical of splice sites 
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This Table is compiled from 169 and 179 examples of 5' and 3' ess 
predictions, respectively, that are supported by only single ESTs 
(Supplementary Table SI). The relative proportions of the bases T, 
C, G and A are shown as a percentage at five positions both 
upstream and downstream of the predicted cryptic cleavage site. The 
most frequently occurring bases are shaded. 



amount of EST data that is available. Alignment of the 
predicted ess indicates that the large majority are correct 
because they conform to splice site consensus sequences 
(Supplementary Table S4). All of these ess were used at 
low frequencies, similar to the HBB example of Figure 2A. 
Figure 4B shows how we record the position of ess; in this 
example, there is a ess in the human RPS5 gene that 
exactly matches the position of introns in Drososphila 
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Figure 4. Comparison of intron and ess positions for a small part of 
the ribosomal protein gene database (23). (A) An alignment of 14 
amino acids of the RPS5 gene that marks the position of introns for 
23 species. (B) The same alignment including two ess positions 
identified by CSF. (C) An alignment of part of RPL7A that illustrates 
the conservation of two ess that also match an intron in 
Chlamydomonas reinhardtii (Cr). * - marks the position of ess that 
match introns; A - marks conserved ess. 




1 FANCA 161 bp 




Figure 3. Experimental confirmation of CSF predictions. Predicted ess are shown by the vertical boxes for the indicated genes. Active ess would be 
expected to generate PCR fragments of the sizes shown in the gene diagrams. PCR products marked with asterisks were sequenced in order to 
confirm the use of the predicted ess (Supplementary Figure S3). Messenger RNA was prepared from the human cell lines K562 (lane 1); HEPG2 
(lanes 2, 3) and primary mesenchymal stem cells (lanes 4-6) and used for RT-PCR with the indicated primers (see Supplementary Data). 
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(Dm) the honeybee (Am) and a fungus (Cc) and a nearby 
ess in the mouse (Mm) that exactly matches the position of 
introns in three plant species (At, Os, Cr). We have 
identified 135 ess to date and of these at least five are 
conserved between species (Figure 4C and 
Supplementary Figure S4). Thirty-three out of 135 ess of 
the ribosomal gene family exactly match the position of 
introns in other species (Supplementary Figure S4). This 
proportion is smaller than the 11/14 match observed for 
the actin gene family (Supplementary Figure S2), but is 
still highly significant (P = 2.2 x 10 -15 , see 'Materials and 
Methods' section). 



DISCUSSION 

The CSF program is designed to identify transcripts that 
are generated through the low level use of ess by normal 
genes. In addition, CSF also identifies a subset of alterna- 
tive splice sites that are similar to ess, but are used at a 
greater frequency. Both types of splice sites are reported to 
be activated in genetic disease as a result of mutations that 
disrupt the function of nearby competitive splice sites 
(4,5). Our analyses show that ess predictions by CSF are 
very reliable and so we would expect, for example, that 
many more of the predicted ess and alternative splice sites 
of Supplementary Table SI will be discovered to be 
associated with genetic disease. The identification of ess 
by CSF is limited by the amount of available transcript 
data but this will improve as further transcript sequences 
become available, particularly with the advent of mRNA 
deep sequencing (24). 

CSF was designed primarily to identify ess in highly 
conserved gene families such as actin in order to 
advance our understanding of intron origin. We found 
that about 25% of the ess within the coding sequence of 
the large family of genes that encode ribosomal proteins 
exactly match the position of introns that are present in 
equivalent genes from other species (Figure 4 and 
Supplementary Figure S4). This compares to a 78% 
match between actin ess and introns, however, there is 
far more phylogenetic data available for the actin gene 
family and so relatively more introns have been dis- 
covered. Consequently, the ess that are recorded in 
Supplementary Figure S4 predict the positions of as yet 
undiscovered introns in the ribosomal gene family. 

The well-known and valuable early and late models of 
intron origin both assume that the splicing machinery 
evolved for the purpose of removing introns that were 
either present in the most ancestral genomes or were 
inserted after the separation of the prokaryotes (25-27). 
At the time the models were proposed, non-intronic 
splicing information was not generally thought to be of 
major relevance and so is not an important feature of 
either model (28). 

However, it is now established that exon junction se- 
quences are older and better conserved than most if not 
all introns and were sites for intron insertions during evo- 
lution (18,29-33). Indeed, a number of intron properties 
such as their phasing with respect to the coding sequence 
(29-30,34) and their location with respect to protein 



structure have now been largely attributed to the 
flanking exon junction sequence (35,36). 

Consequently, our finding that ess often match the 
position of exon junction sequences in gene homologues, 
strongly indicates that ess were targeted by intron inser- 
tions during evolution. Our data also indicates that the 
information that lies outside some introns is not only 
conserved with homologs that lack such introns but is 
also capable of being independently recognized by the 
splicing machinery and can define the position of the 
'missing' introns (13). This is a striking observation 
because although the splicing information that flanks 
introns contributes to their recognition, there is no mech- 
anistic reason for this information to be autonomous 
rather than auxiliary in nature for the purpose of intron 
removal (37). It is unlikely that a ess is recognized from 
the immediate splice-site information alone (8-10), sug- 
gesting that other information such as splicing enhancer 
sequences are also conserved between some gene homo- 
logues independently of intron presence. 

All of the available evidence indicates that the splicing 
information that flanks introns was largely in place prior 
to their insertion (18,29-33). The key question we have 
started to address here is whether such information 
might have been functional. Our data so far does not 
prove but it certainly supports the suggestion that infor- 
mation of this nature could have been used by the splicing 
machinery for splicing purposes prior to the arrival of 
introns, which is a rather different concept to the early 
and late models of intron origin (28). 

The splicing information that flanks introns is perhaps 
similar to the splicing information that enables the mRNA 
of intronless genes to interact with components of the 
splicing machinery for the purpose of mRNA biogenesis, 
nonsense mediated decay or alternative splicing (38-43). 
For genes that have introns there are still many reports of 
intron-independent splicing between exonic splice sites 
(44-53). If intron-independent splicing is an ancestral 
mechanism then it may be far more prevalent across 
species than is currently realized. 



SUPPLEMENTARY DATA 

Supplementary Data are available at NAR Online. 
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